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We calculate the time evolution of fluctuations for rare particles such as e.g. kaons in f AGeV or 
charmonium in 200 AGeV heavy ion collisions. We find that these fluctuations are a very sensitive 
probe of the degree of chemical equilibration reached in these collisions. Furthermore, measuring 
the second factorial moment the size of the initial population can be determined. 



I. INTRODUCTION 

Statistical models have long been used as a tool to de- 
scribe particle production in heavy ion and in high energy 
particle collisions . Recent analysis have shown that 
these models can indeed give a satisfactory description 
of the multiplicities of most hadrons measured in A-A 
collisions for bombarding energies ranging from 1 AGeV 
(SIS) to 160 AGeV (SPS) |,|. Especially at the low 
energies (SIS) and for the rare particles such as kaons, 
the success of the statistical description has been a puz- 
zle and could not be understood within the state-of-the- 
art transport models For instance a recent analy- 
sis for chemical equilibration within a transport model 
gives chemical equilibration time of the order of 300 
fm/c for a situation relevant to 1 AGeV heavy ion col- 
lisions, considerably larger than the typical duration of 
such reaction, which is of the order of 30-50 fm/c. On 
the other hand the success of the statistical model cannot 
be disputed. Not only are the particle ratios reproduced 
for central collision over a wide range of bombarding en- 
ergy, but also the centrality dependence at low energies 
is consistent with the predictions of the statistical mod- 
els. In particular the almost quadratic dependence of the 
kaon multiplicity on the number of participating nucleon 
at SIS follows directly from the statistical model, once 
strangeness conservation is taken into account exactly . 
Transport models on the other hand fail to reproduce the 
observed centrality dependence and particularly energy 
dependence of the K/ 7r-ratio [|j . 

In a recent paper |^ some of us have shown, that the 
chemical equilibration time is considerably shortened if 
the strangeness conservation is taken into account ex- 
plicitly. However, using the new rate equations derived 
in this work combined with the cross sections as given 
in p4| one still arrives at equilibration time substantially 



exceeding the lifetime of the system. 

Does this mean that the success of the statistical model 
is a pure coincidence? This is very unlikely as statistical 
model naturally explains most of the basic features of 
experimental data in a very broad energy range from SIS 
up to SPS. It is conceivable that there are additional 
processes at work, like e.g. many particle collisions ]l5[ | or 
in medium modifications of hadron properties [ p6[ , which 
are not yet taken into account in transport models. 

Thus, a direct experimental determination of the rate 
of equilibration in heavy ion collisions is called for as it 
would possibly provide evidence for new physics. In this 
paper we will demonstrate that the fluctuations of rare 
particles is a very sensitive probe of the degree of equili- 
bration reached in these collisions. Such a measurement, 
though certainly difficult, could for the first time provide 
a direct experimental evidence for chemical equilibration 
in heavy ion reactions. 

This paper is organized as follows. In the following sec- 
tion we set up the formalism. Then we present the results 
for the time dependence of the second factorial moment 
for several initial conditions. Before we turn to observa- 
tional issues we also will discuss the case when no con- 
straints from any U{1) charge conservation are present. 

II. FORMALISM 

In Ref. 1^ the rate equation for particles which are sub- 
ject to an explicit U{1) "charge" conservation has been 
derived. Considering a binary process 0102 — > 6162 with 
a ^ b one arrives at the following master equation for the 
probabilities P„ to find n particles 
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- [n2p„ - (n + ifPn+i] , 
where n = 0, 1, 2, 3, • • •. Here 

e = G{Na,){Na,)/L, 

and the dimensionless time variable t is defined as 

L 



(1) 



(2) 



(3) 



so that T is measured in units of the relaxation time 
Tq = V/L lH]. The momentum-averaged cross sections 
for the gain process aia2 6162 and the loss process 
bib2 aia2 are defined as G = {(Jgv) and L = (ctlu), re- 
spectively. The ratio of these momentum averaged cross 
sections is related to the ratio of equilibrium particle den- 
sities involved 



L da^a1^K2{aa^)da^al^K2{aa^)'' 



(4) 



where dkS denote the degeneracy factors, and ak = 
ruk/T. 

Eq. (|l|) has no obvious solution and needs to be solved 
numerically. The asymptotic (equilibrium) probability 
distribution, on the other hand, has been derived in |9[] 



P. 



leading to 



g2 ^3 

^^-2 + 3 



\ I cq. 



(5) 

(6) 
(7) 



The above general rate equation is valid for arbitrary 
values of {N) for particle production constrained by C/(l) 
charge conservation. It reduces to the grand canonical re- 
sults for large {N) and to the canonical results for small 
{N). It provides a generalization of the standard rate 
equation beyond the grand canonical limit. It was shown 

that for rare particle production the equilibrium abun- 
dance is much smaller and the relaxation time is much 
shorter than expected from the standard rate equation. 
In this paper we will discuss further consequences of the 
generalized rate equation and in particular study the time 
evolution of the multiplicity fluctuations. We want to 
demonstrate that the combined information on both {N) 
and {N"^) can help to determine the degree of chemical 
equilibration. 



III. RESULTS 



In this work we will be mostly concerned with the fluc- 
tuations of the particle number in the case of rare parti- 
cle, i.e. {N)^^ <C 1 or, equivalently, e <C 1. In particular 



we will investigate the behavior of the second factorial 
moment F2 



^ {N{N-l)) 



(8) 



From Eqs. and (0) in the limit of small e the equilib- 
rium value for F2 is given by 



F2 = i + - 
2 6 



{N) 



cq. 
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(9) 



In order for the second factorial moment to be a sen- 
sitive probe of the degree of equilibrium achieved, one 
needs to investigate its initial value. Here we consider 
two distinct cases: 

1. The initial particle number is considerably smaller 
than the equilibrium value. This is relevant, for 
example, for kaon production in 1 AGeV heavy ion 
collisions. 

2. The initial particle number is considerably larger 
than the equilibrium value. This might be relevant 
for charm production in 200 AGeV heavy ion colli- 
sions 

In the first case, where the initial particle number is 
small, let us consider two scenarios. On one hand, let us 
assume that initially the probabilities P„ are distributed 
according to a Poisson distribution: 



Pn{T 



= -^e 



No 



(10) 



where A^o is the initial average number of particles. In 
this case, the factorial moment obviously starts out at 



F2{t = 0) = 1 



(11) 



and decreases by a factor of two until equilibrium is 
reached. On the other hand one may assume that ini- 
tially there is at most one particle in a given event. In 
this case the initial conditions are 



(12) 



which we will refer to as 'binomial' initial conditions. As 
shown in Appendix^, F2 starts out at F2 = 0, but almost 
immediately reaches a maximum after a time of the order 
of 



Po(r 


= 0) 


= 1-^0, 


Pi(t 


= 0) 


= A^o, 


Pn{r 


= 0) 


= n > 1 



(13) 



and for Nq/N^^, < 1, P^"''^ ~ 1. 

Therefore F2 approaches equilibrium from above and 
a measurement of F2 > 1/2 will indicate the degree of 
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equilibrium that has been reached in a heavy ion coUi- 
sion. The detailed time and height, where F2 reaches a 
maximum depend, of course, on the input parameters A^o 
and e. The dependence of F^'^^ on the ratio No/Ncq. for 
e = 0.1 is shown in Fig. |^ as the full line. The dashed line 
in Fig. |l| is obtained by assuming Tmax = SiVo/A^^eq. show- 
ing that indeed the time scale for reaching the maximum 
is given by No/N^q.. 

We further see that for small A^o the factorial moment 
essentially immediately reaches a value close to = 1, 
giving a factor of two sensitivity on the degree of non- 
equilibrium established in the collisions. Obviously, in 
the case where A'o — e = N^q. this sensitivity is lost, as 
the equilibrium value is very close to the initial value. 



N„/<N>„ =0.01 




FIG. 1. Location and value of the maximum of F2{t) 
for a range of values of No/Ncq. as indicated in the figure. 
Here e = 0.1 has been used. The dashed line assumes that 

Tmax = 3No/Neq. 
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FIG. 2. Time evolution of the factorial moment F2 for sev- 
eral initial particle numbers A^o (thick lines) . The thin lines 
show the result of the approximate formula ( ) . Here e = 0. 1 
has been used. 

Assuming binomial initial conditions in Fig. 0we show 



the full time evolution for several initial particle numbers. 
For small times the approximate solution ( |A2| ) is also 
shown. Clearly, the equilibrium value of F2 is reached 
from above, but the effect becomes small as the initial 
particle number becomes comparable with the equilib- 
rium value. 

In Fig. H we show the time evolution of F2 for differ- 
ent choices of e or, equivalently, for different equilibrium 
particle numbers. Obviously, the larger the equilibrium 
particle number is the smaller is the effect. 
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FIG. 3. Time evolution of the factorial moment F2 for sev- 
eral equilibrium numbers or equivalently values of e. Here 
binomial initial conditions with No/Naq. — 0.1 have been cho- 
sen. 

For the second case when the initial population is much 
bigger than the equilibrium population, then the annihi- 
lation process must dominate during the early stage of 
evolution for a short period of time. In this case, one 
can look for a perturbative solution around e = 0. Since 
there can be many different initial conditions with a large 
initial population, it is better here to use the generating 
function [01 



Kr,a;) = ^x"P„(r), 



n=0 



and the equation it satisfies 

dg 



^ = (1 - x) {xg" + g'- eg) . 



(14) 



(15) 



where the prime indicates a derivative with respect to x. 
The averages needed to calculate the first two factorial 
moments are given by 

(A^) = .g'(r, 1) and {N{N - 1)) = .g"(r, 1). (16) 

Details of the perturbative procedure can be found in 
Appendix |^. To understand the qualitative picture, first 
consider the initial time t — 0. Since we assume (A^) ^ 1 
initially, we must have {N{N—1)) = 0{{N)'^) and hence, 
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F2=C'(l)atT = 0. For r > 1, the first order perturba- 
tive solutions are 



{N) = {N)^ 



2e 



(NiN-l)) = (iV(iV-l)),q. + - |ai|e-- + ©(e"^-), (18) 



where ai is a 0(1) constant determined by the initial 
condition. From the above equations and using Eqs. (|^) 
and (|^) the second factorial moment is then given by 



l £^ + (4/5)e|fli|e-- 
^ (e + |ai| e ^) 



(19) 



For times 1 ^ t ^ — In e, the exponential terms in 
Eq. (|ig| ) dominate. As a result, F2(r) = 0(e) within 
the interval 1 < r < — Ine. Note that for arbitrarily 
small e, this interval can be arbitrarily long. Also, since 
0(e) < i^2(0) and 0(e) < i^l'^', must reach a mini- 
mum somewhere inside that interval. 

For illustration, let us choose a Poissonian initial distri- 
bution with Nq — 5 and also set e = 0.1. The numerical 
solutions and Eq. (|l^) as well as the second order per- 
turbative solution are displayed in Fig. ^ The numerical 
solution clearly shows the rapid initial decrease and the 
subsequent slower rise to the equilibrium value. To illus- 
trate the duration of the small F2 interval, we also show 
the full numerical result with e = 0.001. The longer pe- 
riod of time where F2 stays to be 0{e) is clearly visible. 



Full solution 
Second order solution 
First order solution 
Full solution : E = 0.001 




FIG. 4. The second factorial moment as a function of time. 
Initial distribution is a Poisson distribution with A^o ~ 5. 
e = 0.1. The solid line represents the numerical solution. 
The dashed line is the result of the second order perturbative 
calculation and the dot-dashed line is the result of the first or- 
der perturbative calculation. Also shown is the full numerical 
solution with the same initial condition and e — 0.001. 

To sum up, F2 as a function of time must: 
(i) Start from 0(1). 



(ii) Reach the minimum value of 0(e) at r 

l/2 + 0(e). 



1 which 



is much smaller than ^2'^' 



(17) (iii) Stay 0{e) until r ^ - Ine. 



(iv) Approach the equilibrium value from below after 
r —Ine. This is in contrast to the first case 
where we considered A^o ^ (-^)oq • 

Hence, if experimental value of F2 is smaller than 1/2, 
then it is a strong indication that the equilibrium is not 
reached and furthermore it also indicates that the initial 
population was much larger than the equilibrium one. 



IV. ABSENCE OF ;7(1)-CHARGE 
CONSERVATION 



It is interesting to also study particle production with- 
out the constraint of U{1) charge conservation. Obvi- 
ously the time evolution equation for the multiplicity 
distribution should be different from what we have de- 
scribed so far for strongly correlated processes. Let us 
consider a general process a + b c + d without any 
constraint of charge conservation. We denote P„ as the 
multiplicity distribution for particle c. It then satisfies 
the following evolution equation. 



P„] - [nP„ - (n -f l)P„+i] 



(20) 



where f = t/ (Nd) = tL/V/ (Nd) is a scaled time and 
i = e/ (Nd) = G{Na){Nb)/L/{Nd). 



The generating function g(f,x) for P„ 



ff(f,x) = ^x"P„(f) 

n=0 

satisfies the following partial differential equation, 



df 



{l-x)[g'-eg] 



(21) 



(22) 



It is interesting to note that the above equation is very 
similar to Eq. (15) for a constrained system except that 
it does not contain the second derivative on the right- 
hand side. Therefore, for certain period during which 
one can neglect the second derivative of the generating 
functions, the evolution of the multiplicity distribution in 
a canonical system should be similar to that of a grand 
canonical. 

The general solution to the Eq. (p2|) can be found 



(23) 



with the initial condition g(0,x) = go(l — x). The nor- 
malization condition g{f, a; = 1) = ^ P„ = 1 also implies 
5o(0) = 1. One can readily find the equilibrium solution 
in the limit f = 00, 
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(24) 



with the corresponding equilibrium multiplicity distribu- 
tion 



71, cq. " I , 



(25) 



which is a Poisson distribution with averaged multiplicity 
(iV)cq. = e. One can also easily calculate the first and 
second factorial moments of the multiplicity distribution, 

(N) =e+{{N)o-e)e-' (26) 
{N{N - I)) = {N{N ~ l))oe-^^ 

- 2€{N)oe-^e-^ - 1) + e'{e-^ - if. (27) 

They both approach to their equilibrium values exponen- 
tially. 

One interesting case of a special initial condition is 
go — 1 when the initial multiplicity {N)q is zero. Com- 
paring Eqs.(p3[) and (p^), one finds that the multiplicity 
distribution in this case remains a Poissonian, and con- 
sequently the factorial moment F2 = 1 at all times. For 
the binomial initial conditions (Eq.(p^)), 



(28) 



the factorial moment F2 starts out from F2{f = 0) = 0, 
but approaches the equilibrium value via 

F2^l' [Afo/(eV - 1) + No)]'' (29) 

at a time scale of 

f ~ ln(l + V2No/e) ~ v^iVo/?. (30) 

Comparing the results presented here with the previous 
section it is clear that an additional constraints imposed 
by the U(l) charge conservation laws are implying a cru- 
cial modification of not only the equilibrium probability 
distributions of particle number but also their fluctua- 
tions and time evolution towards the equilibrium limit. 

V. TOWARDS EXPERIMENTAL OBSERVABLES 

So far we have discussed the production of a single 
species of particles with conserved quantum numbers, 
such as e.g. K~^-K~-pa,irs. In reality, however, one has to 
deal with more than one species. For example at 1 AGeV 
heavy-ion collisions the relevant strange degrees of free- 
dom are and if" which carry positive strangeness 
and the A and S hyperons carrying negative strangeness. 
The anti-kaons are not relevant in this case as the ratio 
of K-/K+ is very small (~ 2% at 1.5 AGeV ). The 
results shown in the previous section thus apply to the 
combined multiplicities for and i.e. 



N = Nk+ + N, 



(31) 



Very often, experiments can either measure K'^ or 
but not both species at the same time. The above mas- 
ter equation has been extended for more than one particle 
species carrying the conserved quantum number in ref- 
erence ]lT| . The equation for the combined probability 
Pi,j to find i and j mesons is given by 



dr 



^,3)+<i2^{P^,j-l - Pi,,) 
^1 



^ ei{P,-i^j ~ P, 
-{l{l+j)P,^j - (z + + l)P.+ i,,) 

U + l){i+J + l)P;,j + l), 



1^1 



where r = and ei,2 = Gi^2{Na^){Na^) / Li^2- 
equilibrium solution is given by |0| 



i+3 
'^tot 



poq. _ 



{i+j)\^^ 



(32) 
The 

(33) 



with etot = ci + £2- Note that the equilibrium probabil- 
ity distribution is the product of the distribution of pairs 
according to Eq. (|^) and a binomial distribution, which 
determines the relative weight of the individual particles, 
in our case the and . For the equilibrium config- 
uration the relevant expectation values are then easily 
computed 



(34) 



(^1 )eq. = + - h) (iV>,„, , (35) 



Fl 



equil 



f-tot 



^ ' 2 6 



(36) 



where /i = ei/ctot- The average ^ denotes the aver- 
ages given m Eqs. (|) and (|) with e — > Ctot For small 
ei_2 the effect of the second species only appears at next 
to leading order in i^2- 

The master equation governing the particle of interest, 
say the K'^ is obtained by summing equation (32) over 
the index of the particle which is not observed. 



dr 



e^{P,^,,-P,) 

-{ep,-{i + i)'p,+i) 



with 



P,. 



(37) 



(38) 



Comparing with the original equation (^, the presence 
of the other species, the in our case, leads to the last 
two terms of equation (|37|). However, in the situation 
of interest here, where N^o ^ 1 these terms can be ne- 
glected. Thus we recover the original equation governing 
the time evolution of the . This can also be seen in 
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Fig. where we compare the evolution of F2 based 
on Eq. (|l|) with that based on Eq. (|3^). For the case 
at hand, namely kaon production in heavy ion collisions, 
isospin symmetry suggests that the production and ab- 
sorption cross-section for K'^ and are roughly the 
same, i.e. Li — L2 and ei — 62. Here we have assumed 
that due to isospin symmetry the collision rates for both 
kaon species are identical. 

To summarize the effect of the second species leads to 
a sub-leading correction ~ cko /6 ~ {K^) /6 which in the 
present context is negligible. 
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FIG. 5. Comparison of time evolution of /2 for one and 
two particle species assuming binomial initial conditions and 
A^o = 0.01 

Let us now turn to the question of how to measure 
these fluctuations in experiment. Since we arc consid- 
ering rare particles, the measurement of fluctuations are 
obviously difficult. Here we propose to study the second 
factorial moment. Therefore, besides the inclusive par- 
ticle number {N) one also has to measure {N{N — 1)). 
The latter expectation value, however, is directly related 
to the two particle density 



P2{Pl,P2) 



(Ppi d^p2 



{N{N - 1)) = / d^pi d^p2 P2{P1,P2). 



(39) 



(40) 



It is interesting to note (see also [Q ) that the same 
information enters the measurement of so called HBT or 
Bose-Einstein (BE) correlations . The BE-correlation 
function as a function of the relative momentum is de- 
fined as 



CsEiq) 



P2{q) 

Pii(9)' 



(41) 



where 

P2{q) = I dpidp2P2{pi,P2)Si\pi -P2I - q) (42) 



and 

piiiq) = J dpidp2Pi{pi)p2{p2)K\pi -P2I - q) (43) 

Usually the Cbe is parameterized as 

CBB(g) = l + Ae-«'«' (44) 

so that outside the correlation region q ^ 1/ Rsource the 
correlation function assumes a value of one. However, 
in case of rare particle subject to U{1) conservation law, 
i.e. kaons at low energy heavy-ion collisions, this will be 
different: 

In terms of p2{q) and pii{q) the factorial moment F2 
is given by 



Idqp2{q) 

fdqpiiiq) 



1 

2' 



(45) 



which, as shown, assumes a value of = 1/2 in equilib- 
rium as a result of strangeness conservation. Using, on 
the other hand, the standard parameterization for the 
BE-correlation function, one obtains 



P2(g) = Pii(g)(i + Ae-*'^') 

and hence 

Jdqp2{q) 



(46) 



Fo 



Jdqpiiiq) 



1 + A / dqe 



p,liq) > 1, (47) 



where the second term is only a very small correction of 
the order of a few percent, given a typical source size of 
Rsource — 5 f m . Obviously the standard parameteri- 
zation for the BE-correlation function is not adequate for 
the case of rare particles subject to a conservation law. 
Since the correlations due to strangeness conservation are 
not expected to introduce any momentum dependence, 
we thus predict, that the BE-correlation for rare parti- 
cles subject to C/(l) conservation should asymptotically 
approach a value of 1/2, i.e. 



Ci3B(<Z» l/Rs 



(48) 



Therefore, the rather difficult measurement of kaon 
correlations at heavy-ion collisions at 1 AGeV would not 
only provide information about the size of the kaon emit- 
ting source but, more importantly, would also be a di- 
rect measurement of the degree of equilibrium reached in 
these collisions. 



VI. CONCLUSION 

In this work, we addressed the question of number fluc- 
tuations of rarely produced particles. We carried out the 
analysis by solving the master equation derived in Ref. 
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IPI . The most important aspect of the master equation 
(|l|) is that we can treat the conservation laws that govern 
the rare particles exactly. For instance, the strangeness 
conservation for kaon production at the SIS energy can 
be treated in this way. Comparing the results of sec- 



tions [II and IV certainly shows the difference such a con- 
straint makes. In previous papers, some of us explored 
the consequence of requiring the exact conservation on 
the behavior of the average multiplicity in equilibrium 
as well as in evolving systems. In this work, we investi- 
gated the time evolution of the second factorial moment 
F2 = {N{N - 1)) / {Nf to explore the possibihty of us- 
ing it as a non-equilibrium measure. 

To cover a wide range of physical phenomena, we stud- 
ied two extreme cases, (i) The initial population of the 
rare particle is much larger than the equilibrium popu- 
lation, (ii) The initial population is much smaller. Our 
main conclusion is that the measurement of F2 can cer- 
tainly tell us if the equilibrium has not been reached. 
Moreover, the approach of the second factorial moment 
towards the equilibrium depends very much on the initial 
condition. Assuming that the equilibrium population e is 
very small, we see that the smaller initial population re- 
sults in the approach from above to the equilibrium value 
(Fj*^' = 1/2 -I- 0{e)). On the other hand a larger initial 
population results in the approach from below with a 
long period of very small i^2- Hence, the experimental 
value of F2 can immediately tell us if the initial pop- 
ulation was smaller {F2^^ > 0.5), larger (FI'^p < 0.5, 
possibly F2^^ ^ 0.5) or the system has already reached 
the equilibrium before the freeze-out (^'2'^'' ~ 0.5). 

In summary: the essential point of this work is that 
in a system in thermal equilibrium the average particle 
number fixes the fluctuations. In case of rare particles 
subject to U{1) charge conservation these fluctuations 
are different from the simple Poisson type finite number 
fluctuations and thus provide a measure of the degree 
of equilibration reached in a system. This measurement 
can either be achieved by measuring the inclusive one and 
two particle densities or via the well known Bosc-Einstcin 
correlations. 



APPENDIX A: INITIAL TIME EVOLUTION FOR 
SMALL INITIAL PARTICLE NUMBERS 

The leading time dependence of the probabilities P„ (r) 
can be obtained by Taylor expansion around the initial 
time r = 



PniT)=Pnir = 0)+J2 



1 d"Pn 



(Al) 



The time derivatives can be obtained by iteratively ap- 
plying Eq.(|l|). To order one obtains 



{NiN - 1)) = e'a' { (2t„ + rl) - a{hTl - -r^) + 0{o?) 



{NY 



2 2 
e a 



((1 + T^f - «(2t„ + 3t2 + t3) -I- 0{a^)) 

(A2) 



where we have neglected higher orders in the small vari- 
able 



< 1. 



We have also rescaled the time Tq. according to 



(A3) 
(A4) 

(A5) 



Initially, at r = 0, the factorial moment starts out at zero 

(A6) 



F2(.^0)^M^^0. 



{NY 

However, after a very short time of the order of r = 
corresponding to Tq, = 1 the factorial moment assumes a 
value 



F2(r. = l)^- (A7) 
which is larger than the final equilibrium result 

F2(T^C50)~i. (A8) 
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Therefore, one expects that the factorial moment F2 
reaches a maximum value of about i^2 — 1 after time 



of the order of r = 



No 



Furthermore, F2 approaches 



equilibrium from above thus a measurement of -F2 > 1/2 
will indicate that equilibrium has not been reached in a 
heavy ion collision. 



APPENDIX B: PERTURBATIVE SOLUTION 



To solve (Eq.(^) 

1^ = (1 - x) (xg" +g' - eg) 



(Bl) 
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perturbatively, we first make an ansatz 

oo 

g{T,x) = g^q {x) + ^ e"" ^a„ft,„(a 

with 



(B2) 



(B3) 



By substituting the expression (B2) into Eq.(Bl) and 
cohecting the same powers of e, one can easily show that 



the equation for ft,!"' is 

(1 - X) {xhiT + h^^^') + n'h^^^ = 
which has the solution 



hi^\x)=F{n,-n;l;x)^J2 



s=0 



is 



1^2 



(B4) 



x' (B5) 



where F{n,—n]l;x) is a hypergeometric function. All 

is) 

other hn s are determined by the functional relations 

- n'hi^^ = {l-x) (xh\^^" + h^' hlr^^) (B6) 

Keeping terms up to e and e~^'^ yields the following 
expressions for the relevant averages 



(N) =g'(r,l) 



Icq. 



aie 



. 2e\ 



At 



(B7) 



{N{N-l)) =g"{T,l) 



(B8) 



To make a statement on how the averages behave as a 
function of time, one needs to know the coefficients a„. 
Obviously, the initial condition fixes these coefficients 

oo oc 

g{0, x) = geq.(x) + ^ a„ hr,{x) = ^ x'Ps{Q) (B9) 

n=l s=0 

In general, to get a„, one needs to invert 

oo 

P,(0) = ^^„a„+0(e) (BIO) 

n—s 

where 



^ / o — n 



2 2n 
I — 71 



(Bll) 



It is not a trivial task in general to solve the above equa- 
tion for all a„. A simple procedure to solve for a„'s can 



be given only if we ignore the C(e) corrections and also 
if there is a last index N for which Pn is non-zero. In 
that case, we can write 



N 



(B12) 



This is a triangular linear system of equations and can 
be easily solved by first getting = Pn{0)/Ann- and 
then QN-i and so on. 

To have a general understanding of how a^'s behave, 
first consider the condition 



oo oo 



s=0 



sn 0,n 

s—0 n—s 

oo n 



(B13) 



ri=0 s=0 



From Eqs. (|B4|), ( |B5| ) and (|B11|), it is easy to see that 

n 

J2Asn^hl°\l)^Son (B14) 
s=0 

Hence, the above condition constrains 

ao = 1 (B15) 
Otherwise, a„'s only have to make < -Ps(O), or 



(B16) 



If s is odd then A^n < for all n > s and if s is even 
then Asn > for all n > s. To keep the probabilities 
positive, a„'s must have alternating signs starting with 
fli < and |a„| must be a monotonic decreasing function 
of n. Furthermore, to keep the probabilities finite, |a„| 
must decrease faster than any power of n. 

Numerical investigation shows that the size of ag is 
0(1) up to s ~ \//Vo and from then on \as\ falls like a 
Gaussian (faster for larger A^'o 15). Empirically, 



an = (-l)"2cxp (-nVMo) (" > 1) 



(B17) 



where Mq « Nq works up to A'o ~ 15. For larger Nq, 
some small P, ca n bec ome negative. With Mq = 15, the 
above expression (B17) gives 



ai = -1.87101 02 = 1.53186 ag = -1.09762 

Qi = 0.688308 as = -0.377751 (B18) 

These coefficients result in a probability distribution with 
(A^)^^o = 14.8 and (A^(A^-l))^^o = 210.2 as shown in 
Fig. I 

Using the initial distribution given by a Poisson distri- 



bution with Nq = 15 and solving Eq. (B12) yield 
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ai 



-1.86667 02 = 1.52 03 = -1.08444, 



04 = 0.682074 as = -0.380919 (B19) 
These of course result in {N{N-1))^^„ = = 225. 
0.2 
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FIG. 6. The filled circles repres ents the probability distri- 
bution resulting horn using Eq. (B17) with Mo — 15. The 
open circles represents a Poisson distribution with A^o = 15. 
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For completeness we also quote the results with a Pois- 
sonian initial distribution with A^o — 5, solving Eq. ( B12| ) 
yields 

ai = -1.6027 02 = 0.871375 03 = -0.347491 
a4 = 0.107955 ag = -0.0272918 (B20) 

These values are used to calculate the perturbative solu- 
tions shown in Fig. ^. The Gaussian formula gives 

fli ^ -1.63746 02 = 0.898658 03 = -0.330598 
04 = 0.0815244 05 = -0.0134759 (B21) 
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